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Abstract. The Householder reduction of a member of the anti-symmetric Gaussian unitary ensemble gives 
an anti-symmetric tridiagonal matrix with all independent elements. The random variables permit the 
introduction of a positive parameter /5, and the eigenvalue probability density function of the corresponding 
random matrices can be computed explicitly, as can the distribution of {qi}, the first components of the 
eigenvectors. Three proofs are given. One involves an inductive construction based on bordering of a 
family of random matrices which are shown to have the same distributions as the anti-symmetric tridiagonal 
matrices. This proof uses the Dixon- Anderson integral from Selberg integral theory. A second proof involves 
the explicit computation of the Jacobian for the change of variables between real anti-symmetric tridiagonal 
matrices, its eigenvalues and {qi}. The third proof maps matrices from the anti-symmetric Gaussian j3- 
ensemble to those realizing particular examples of the Laguerre /3-ensemble. In addition to these proofs, 
we note some simple properties of the shooting eigenvector and associated Priifer phases of the random 
matrices. 



1. Introduction 

Gaussian ensembles of random matrices are best known for their application to quantum mechanics. The 
aim of this application (see e.g. [T3]) is to predict the statistical properties of highly excited energy levels 
when the underlying classical mechanics is chaotic; toward this purpose, the Hamiltonian is modelled by a 
random Hermitian matrix H . Time reversal symmetry requires that the elements of H be real, while there 
being no preferential basis in determining the spectrum requires that the joint probability density function 
(PDF) of H be invariant under conjugation by orthogonal transformations. 

These constraints are all satisfied by the Gaussian orthogonal ensemble (GOE), which consists of real 
random matrices H = ^ (X+X T ) , where Xisannxn Gaussian matrix with all entries independent standard 
normals. The joint distribution of all independent entries is thus seen to be proportional to exp(— Tr H 2 /2), 
which is invariant under conjugation by orthogonal transformations, H i— > OHO T . This key property of the 
GOE explains the adjective "orthogonal" in its name. 

In the theory of quantum conductance through a normal metal - superconductor junction the matrix 
modelling the Hamiltonian must have a block structure to account for both electrons and holes. In the case 
when there is no time reversal symmetry, nor spin-rotation invariance, this block matrix must be of the form 

i 



A_ B 

-B -A 



a = aK b = -b 



T 



Such matrices are equivalent under conjugation to i times a real anti-symmetric matrix, and so motivate the 
consideration of anti-symmetric Gaussian matrices H = |(X — X T ) where X is an N x N real Gaussian 
matrix with entries N[0, 1/2]. Here and throughout the paper N[/x, a 2 ] refers to the normal distribution with 
mean \i and variance a 1 . This class of random matrices — referred to as the anti-symmetric Gaussian unitary 
ensemble — has also received recent attention for its appearance in the study of point processes relating to 
the tiling of the half hexagon by three species of rhombi, and to classical complex Lie algebras [9, 5 . 

It is our objective in this paper to initiate a study of anti-symmetric Gaussian matrices from the viewpoint 
of the underlying tridiagonal matrices. To appreciate the possible scope for such a study, let us recall that 
the GOE has a "sibling" ensemble of tridiagonal matrices which shares the same eigenvalue PDF. The latter 
was constructed by applying a well-known numerical algorithm for eigenvalue computation to the GOE. For 
a general n x n real symmetric matrix, a common preliminary strategy in numerical eigenvalue computation 
is to first conjugate by a sequence of Householder reflection matrices so as to obtain a similar tridiagonal 
matrix. In the case of GOE matrices, it was shown by Trotter [18] that the resulting tridiagonal matrix 

l 



exhibits a remarkable property: all elements are again independently distributed (subject only to symmetry). 
Explicitly, one obtains tridiagonal matrices 



(1.1) 



N[0,1] Xn-l 
Xn-l N[0, 1] 
Xn-2 



Xn-2 
N[0, 1] 



Xn-3 



X2 N[0,1] xi 
Xi N[0, 1] 



where N[0, 1] refer to the standard normal distribution and x.k denotes the square root of the gamma distri- 
bution r[fc/2, 1], the latter being specified by the p.d.f. (l/r(fc/2))u fc /2-i e -« ; u > 0. 

Two other Gaussian ensembles also apply to the highly excited states of chaotic quantum systems, namely 
the Gaussian unitary ensemble (GUE) of complex Hermitian matrices and the Gaussian symplectic ensemble 
(GSE) of Hermitian matrices with real quaternion entries (see e.g. [5]). It was pointed out by Dumitriu and 
Edelman [7] that applying the Householder transformation to these matrices gives the tridiagonal matrix 
(jl.ip but with the replacements 



(1.2) 



Xk >-> Xpk 



(k = l,...,n- 1) 



where (3 = 2 for the GUE and (3 = 4 for the GSE. Moreover, with the replacements (|1.2p for general (3 > 
it was proved in [7] that the eigenvalue PDF of the corresponding random tridiagonal matrices is equal to 



(1.3) 



where 



n 

l<j<k<n 



G p , n = (2^)™/ 2 Yl 



r(i + (j + 1)13/2) 
r(i + /3/2) ' 



Our specific aim then is to give random tridiagonal matrices whose eigenvalue PDF generalizes that of the 
anti-symmetric Gaussian matrices (as well as the ensemble consisting of pure imaginary elements, another 
ensemble that can be constructed out of pure imaginary real quaternion elements). As in we find 

that the tridiagonal matrices in question have independent elements (up to the anti-symmetry condition), 
and that these elements allow for a (3 generalization. Explicitly, our study focuses on the family of random 
anti-symmetric tridiagonal matrices 



(1.4) 



X(n-l)0/2 
-X(n-l)0/2 X(n-2)P/2 

— X(n-2)/3/2 X(n-3)/3/2 





'X/3/2 



X/3/2 




For n even, the eigenvalues of these matrices come in pairs {H^j}j=i n/2i > 0j while for n odd zero is 
a simple eigenvalue, and the remaining eigenvalues come in pairs {±iA_y}j = x,...,( n -i)/2! Aj > 0. 

Our main results are given in Theorem 11.11 and Theorem 11.21 for which we give three proofs; each proof 
uses different tools and highlights different properties of the matrices in question. 

Theorem 1.1. With the notation above, the PDF of the positive variables {Xj}, ordered non-decreasingly 
(Xx < A2 < . . • < X n ), for n even is given by 



/2 



(1.5) 



/3/2-l e -A? 



n ( a , 2 -^)^ ; 

l<j<fc<n/2 



while for n odd it is given by 

(i.6) ^- (n nV 2 -v^ n 

13,11 i=l l<j<k<(n-l)/2 

The normalization constants are given by 

n/2 



13 



(1.7) 



Yj r(2j/3/4)r((2j - l)/3/4) 

^." = 11 2T(^ ' 

C/3,„ = r j"// 4 p (»-l),/3i ? ^ odd 



These PDFs are related to /3-Laguerre ensembles (see [4], [7]). More precisely, using the notation of 
[7] , the squares of the eigenvalues of antisymmetric Gaussian /3-ensembles of size n are distributed like the 
eigenvalues of /3-Laguerre ensembles of size [§J and parameters a — for n even, and a — 2fi for n 

odd. 

In addition to the PDFs, we give the distributions of the first components of the eigenvectors. 

Theorem 1.2. Let qj, j = 1, . . . , [§J be the (positive) first component of the independent eigenvector 
corresponding to {Aj}. In addition, if n is odd, let z be the (positive) first component of the eigenvector 
corresponding to the (simple) eigenvalue. Then 

• if n is even, the vector (qi, . . . , q n /2) has Dirichlet distribution Z? n / 2 [(/3/2)™^ 2 ] (here (/3/2)™/ 2 denotes 
(3/2 repeated n/2 times); 

• if n is odd, the vector (qi, . . . , qi n -x) /2> z ) has Dirichlet distribution D^ n _iy 2 [((3/2Y n ^ 1 ^ 2 , [3/4]. 

The rest of the paper is organized as follows. In Section [2] we begin by recalling the definition of the 
two classical ensembles of anti-symmetric Gaussian matrices, and show that the Householder reduction of 
the anti-symmetric Gaussian unitary ensemble gives (|1.4|) with (3 — 2. Although we do not give details, the 
Householder reduction in the case of pure imaginary real quaternion elements gives (jl.4|) with (3 = 4. 

We show three different approaches to the eigenvalue PDF computation for the matrix model (|1.4|) with 
general (3; all of them extend to the computation of the PDF of the first component of the corresponding 
eigenvectors. Each proof is self-contained. The first one, in Sections [3] and |4j uses the characteristic poly- 
nomial of the matrix, which is constructed inductively. A second approach, shown in Section [5l is based on 
Jacobians of tridiagonal anti-symmetric matrix. A third proof is presented in Section [6l where we compute 
the PDF of Theorem 11.11 by use of an orthogonal similarity transformation between the anti-symmetric [3- 
ensemble and a /3-Laguerre ensemble. Finally, in Section [7] we make note of some simple properties of the 
shooting eigenvector and associated Priifer phases of the random matrices (|1.4I) . 

2. Householder reduction 

As mentioned, there are two classical ensembles of anti-symmetric Gaussian matrices. One is the anti- 
symmetric Gaussian unitary ensemble, which is specified in the paragraph below (|1.3| . The other is the 
anti-symmetric of interest in statistical physics might have /3-generalizations; one case in point is that of the 
antisymmetric Gaussian of interest in statistical physics might have /3-generalizations; one case in point is 
that of the antisymmetric Gaussian of interest in statistical physics might have /3-generalizations; one case 
in point is that of the antisymmetric Gaussian of interest in statistical physics might have /3-generalizations; 
one case in point is that of the antisymmetric Gaussian of interest in statistical physics might have (3- 
generalizations; one case in point is that of the antisymmetric Gaussian Gaussian symplectic ensemble. A 
member A of this ensemble is formed out of an n x n Gaussian random matrix Y with quaternion real entries 



Y 
1 'j 



z w 
—w z 



(as a complex matrix, Y is 2n x 2n), according to A = \(Y — Y T ). Both the z and w variables are 
required to be purely imaginary, chosen from «N[0,1]. This ensemble has the special feature that the eigen- 
values are doubly degenerate. Generally the eigenvalues of anti-symmetric Hermitian matrices come in pairs 
{±i\j}j = i^_ i [ n /2] (Aj > 0) and with zero as a simple eigenvalue for n odd. 
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The PDFs of the positive eigenvalues for both the anti-symmetric GUE and GSE are known (see e.g. [5]). 
For the anti-symmetric GUE the PDF is proportional to 



Afc) 2 , 



(2.1) 



n/2 

II 

i=l 

(n-l)/2 

n ^ n ^ 

i=l l<j<k<(n-l)/2 

while for the anti-symmetric GSE it is proportional to 



n 

l<j<k<n/2 



2\2 



n even 



n odd 



(2.2) 



n/2 

i=i 

(n-l)/2 

n 

»=i 



n (a. 

l<j<fe<ri/2 



Afc) 4 , 



n even 



Afe" 



n 

l<3<A<(n-l)/2 



Here we give the explicit form of the Householder reduction of anti-symmetric GUE matrices, and so deduce 
a random tridiagonal matrix with eigenvalue PDF proportional to (|2.ip . 

As for real symmetric matrices |20j . in the case of anti-symmetric matrices [fly]ij = i r .., n , the House- 
holder tridiagonalization consists of constructing a sequence of symmetric real orthogonal matrices 
(j = 1, . . . , n — 2) such that 

(n-2) (n-3) . . . q(1) A q(1) q(2) . , , q(u-2) _ S (n-2) 

where _B(™~ 2 ) is an anti-symmetric tridiagonal matrix. 

This process is based on the fact that any column vector x — (xi, . . . , x n ) T can be mapped into ||x||2ei, 
with ei being the first column of the identity matrix I„, by the symmetric orthogonal transformation H = 



'IMP 



where v — x + X\e\ and I \v\ 



\v\\2- The matrix H is called the Householder reflector for x. 



For our antisymmetric matrices, the matrix is chosen so that :— O^AO^ is tridiagonal with 
respect to the first row and first column; the matrix is chosen so that :— O^B^O^ is tridiagonal 
with respect to the first two rows and first two columns, etc.. For this must be of the form 



(2.3) 



U) 



R(n-j) 



where Ry 1 ^ is the Householder reflector corresponding to the j + 1st through nth entries in the jth column 



of oCj-1) . . . OWAOM ■ ■ ■ O^'- 1 ) 
It follows that 



in 



Here 



an, 



"IXJ 



is the i x j matrix of all zero entries. 



( J B«) 1 2=-(B (1) ) 21 = (<*12 + 



1/2 



(2.4) 
etc.. 

With real anti-symmetric matrices formed by — i times an anti-symmetric GUE matrix, it follows from 
(|2.4p that (B^)i2 has distribution equal to the square root of T[(n — l)/2, 1]. 

Anti-symmetric GUE matrices are distributionally invariant under conjugation by an (independent) or- 
thogonal matrix; this is an immediate consequence of the fact that a column vector of i.i.d. Gaussians does 
not change its distribution when left multiplied by an (independent) orthogonal matrix. 

This invariance under conjugation, as well as the structure (|2.3j) of 0^ x \ tell us that the sub-block of B^> 
obtained by deleting the first row and first column is —i times an anti-symmetric GUE matrix of size n — 1. 
Proceeding inductively and remembering the structure (|2.3p . we see that (|1.4I) holds in the case = 2, 



3. Method I part (i): bordered matrices and an inductive construction 

An inductive construction of the tridiagonal matrices (|1.1[) . (|1.2p involving the operations of bordering 
was given in [10] . In [9] this same bordering procedure was used to compute the joint eigenvalue PDF of an 
anti-symmetric GUE matrix and its successive principal minors. Here the working of |9J will be extended to 
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provide an inductive construction of a family of random matrices with eigenvalue PDFs given by Theorem 
11.11 Additional workings from [10] will then be adapted to deduce the tridiagonal matrix Ijl.ip itself. 



Proposition 3.1. Let w be an n-component real column vector of the form 
(3.1) w = 



(wi,wi, . . . ,w n / 2 ,w n / 2 ), n even 

(wi,u>i, . . . ,i0( n _i)/2, W( n -i)/2,b), n odd 

Let 2wj have distribution T[/3/2, 1] and let b 2 have distribution T[f3/A, 1]. Define the sequence of Hermitian 
random matrices A\ = 0, Ai, A3, ... , where each A^ is k x k, by the inductive formula 



(3.2) A 



77+1 



diag A n iw 
-iw T 



where diag A n is the diagonal matrix formed from the eigenvalues of A n . The eigenvalues of each Af. come 
in plus/minus pairs {±Xj}j—i ,{k/2]> Xj > 0, and for k odd zero is also a simple eigenvalue. (Such pairing 
is to be taken as implicit in diag A n , and when n is odd the zero eigenvalue is to be listed last.) Furthermore, 
with the characteristic polynomial of A n denoted by P n (x), for n even we have 

(3.3) ^M.l-V^, 
while for n odd 

P n+ i(x) _ b 2 ( "^ /2 2wf 
[ ' ' xPJx) x 2 ^ x 2 -X 2 ' 

Proof. We see that if v is an eigenvector of (j3.2p with eigenvalue A, then v is an eigenvector with eigenvalue 
—A. Thus, as claimed, the eigenvalues come in plus/minus pairs, and for k odd there is a zero eigenvalue 
corresponding to an eigenvector with all components real. To establish the equations for the characteristic 
polynomial, we first note that by induction we must have 

(Ai, -Ai, . . . , A„/ 2 , -A n/2 ), n even 

(Ai,-Ai,...,A(„_i)/ 2 ,-A( n _i)/ 2 ,0), n odd 

Recalling (|3.ip . it follows from (|3.2p that for n even 

77/2 ^ 

det(xl n+ i - A n+ i) = det(xl„ - A n )(x - ^^w 2 ^- — 



(3.5) diag^„ = 



while for n odd 



det(xl n+ i — A n+ i) = det(icl n — A n )(x — — - ^ tr 



Xj x + A,- 
j=i J J 



2 (7i-l)/2 
\ ^ 2 



x ^— ; ■ , \x — Xj ' x + A, 
These are the equations (|3.3p and (|3.4|) respectively. □ 



The equations (|3.3p . (|3.4|) can be used to compute the conditional eigenvalue PDF for A n+ i, given the 
eigenvalues of A n . We see that the task is to compute the density of the zeros of the random rational 
functions therein. This can be done by appealing to the following known result. 

Proposition 3.2. |101 Corollary 3] Consider the random rational function 
(3.6) R(x) = 1 - V 

7—1 

where each qi has distribution T[si, 1]. This function has exactly n roots, each of which is real, and for given 
{a^ these roots have PDF 



supported on 

X\ > ai > ■ ■ ■ > x„ > a n . 

In the case of n even, the RHS of (|3.3p corresponds to (|3.6p with n i— » n/2, x i— > x 2 , a, i— > A| and each 
distributed as T[/3/2, 1]. In the case n odd, the RHS of (|3.4p corresponds to (|3.6p with n i— > (n + l)/2, 
.x i— > a; 2 , ai i— > Af (i = 1, . . . , (n — l)/2), a(„+i)/2 = 0, qi distributed as r[/3/2, 1] (i = 1, . . . , (n — l)/2) and 
Q(n+i)/2 distributed as T[/3/4, 1]. The sought conditional PDFs can therefore be made explicit. 

Proposition 3.3. Forn even, the PDF of the positive eigenvalues of A n+ i, given that the positive eigenvalues 
of A n are Ai, . . . , A„/ 2 ; * s equal to 



2 n/2 JJ 



(37) tl TT (Xt ~ Xj) TT \x 2 A 2 ^/ 2 - 1 

[ ' (r(/3/2))»/ 2 H ( A 2_ A 2 )/3 -i 11 ^ J 1 

where 

(3.8) xx > Ai > • •• > ic n /2 > A„ /2 . 

for n odd, the PDF of the positive eigenvalues of A n+ \, given that the positive eigenvalues of A n are 
Ai, . . . , A(„_x)/2, is equal to 

(n-l)/2 



(r(/3/2))(»-D/ 2 r(/?/4) 

(n+l)/2 

II X*l' 11 W ^jV (n+l)/2 (n-l)/2 

(3.9) X— !=± l<i<i<(n+D/2 TT TT | 2_ ^0/2-1 

(«-i)/2 n a 2 - a 2 v j -! 11 11 

11 l A i7' l<»<J<(n-l)/2 
1=1 

w/iere 

(3.10) xi > Ai > • • • > X( n _i)/ 2 > A(„_x)/ 2 > X(„ +1 )/ 2 > 0. 
Let the conditional PDFs of the above proposition be denoted 

Gn+l({Zj}j=l,...,[(n+l)/2y,{^j}j=l,...,[n/2]), 

and let the marginal PDF of the positive eigenvalues of A n be denoted Pn({ x j}j=i [n/2])- Then we must 
have that 

Pn+l({Xj}j=l,...,[(n+l)/2]) 

(3-11) = / dXi--- dX[ n/2 ] Gn+l({a;j} 3 -=l,...,[(n+l)/2];{Aj} :! -=i,...,[„/2])Pn({Aj}j=i [ n /2])- 

Jo Jo 
Furthermore, from the definition of Ai we have that 

(3.i2) P2{x) = _L_ x m-i e -*\ x>0 

so p. lip uniquely specifies {p n ({xj}j=i ,...Jn/2])}«=3,4,...- As such, it can be used to verify that the explicit 
functional form of p n is given by Theorem 11.11 

First proof of Theorem 11.11 We make a trial functional form 

, [n/2] 

Pn{{\ i } j =x,..., W 2 ] ) = ^- n e_A * ^ n ( a ? - a ')" 

i=l l<j<fc<[«/2] 
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where C n ^ and a = a n ^ are to be determined. Substituting in (|3.11[) then gives, for n even, 

-^(r(/3/2))«/ 2 YlxT^' 1 n {x] - xlf-' 

°' i+1 '' 3 i=l l<j<fc<n/2 

n/2 n/2 n/2 

(3.i3) = 2 «/ 2 / rfAx-.-dA^ n^n^-*?) n w-^ 8 - 1 

•^"/a fc=l i<j i,j=l 

where R n /2 is the region (|3.8j) . while for n odd 

(n+l)/2 

(r(/3/2))("- 1 )/ 2 r(/3/4) ff x ^+^-^ 

2C„+i !/3 ax 



(n+l)/2 



X 



JJ ^,,+1-/3/2 -Q ( x 2_ 4)/3 -l = 2 (n-l)/2 /■ dAl ... dA( „_ 



(n-l)/2 



i=l l<j<fe<(n+l)/2 x >>-l)/2 



(n-l)/2 (n-l)/2 (n+l)/2 (n- 1)/2 

(3.i4) x n A r- +2 - 3/3/2 n ^-^ n n m-w- 1 

l—k i<j i—l j — 1 

where R[ n _iy 2 denotes the region of integration (|3.1Qp . 

To evaluate the integrals, we make use of the Dixon- Anderson integral (see e.g. [SJ Ch. 3]) 

/n n+1 
dXf-dXn J] (A, A, • J J II A 

- r l<j<fc<n 
n+1 

n r(«o 

where X is the domain of integration 

ax > Ai > d2 > Aa > • • • > A n > a n+ i. 
After the simple change of variables Xj i— > A 2 and the replacements a, i— » x 2 , we see that for 
(3.16) a n>p = 3/3/2 - 1 (n odd) 

the RHS of ([3~T4]) is equal to 

r((n + l)/?/4) 11 ^ ^ ' 
The fact that the above is identical to the LHS of (|3.14[) provides 

(3.17, c,. s gaaaffl ., (nodd) 

and 

(3.18) Qfn+1,/3 = /?/2 — 1 (n odd) . 

Furthermore, (|3.15p with a n+ i = allows us to compute the RHS of (|3.13[) as equal to 



r(K i/3 + l)/2)(r(/3/2))«/ 2 Jr^-i+a.,, tt ( 



~ x k) 



r(n/3/4+K !/3 + l)/2) 11 • L£ > 

This is identical to the LHS provided that 

(3.19, ^ ^m±^m = 1 ( „ev») 

k n +i,/3 1 ((a«,/3 + l)/2) 



and 

(3.20) a n+i,p = ct n ,fi + (3 (n even). 

We observe that (|3.16p . (|3.18[) and (|3.20|) overdetermine a n ,p, as p.!6[) and (|3 . 18[) determine a n ,f3 in all 
cases. Substituting (|3.16p into the LHS and (|3.18p into the RHS of (|3.20[) shows this latter equation to be 
identically satisfied. Further, we observe that (|3.18[) is consistent with the requirement of the initial condition 
(|3.12|) . With otn t p thus specified, substituting in (|3 . 19f> simplifies that formula to read 

,„ n Cn£ r((n + l)/3/4) 

(3-21) - — ttj- = 1 (n even). 

Cn+1,/3 1 (P/4) 

The equations (|3.17|) . (|3.21| together with the requirement C2.fi = r(/3/4)/2 as read off from (|3.12p give 
(|1.7|) for the normalizations. □ 

The values of the normalizations (|1.7[) are subject to an independent check. Thus a well known corollary 
of the Selberg integral (see e.g. [HQj]) gives 



/>oc />oo 71 

W a ,f), n ■= / dx!--- dx n Y[ x i e ~ Xi II \%k-Zj 

J° ^° i=l l<j<fc<n 



If 



n-1 

n 

3=0 



r(l + (j + l)/3/2)r(a+l+j/?/2) 

r(i + /3/2) 



On the other hand, a simple change of variables in the definitions of Cp^ as implied by Theorem 11.11 
remembering too that the eigenvalues therein are assumed ordered, shows 

This is consistent with (|1.7p . 

4. Method I part (ii): A random corank 1 projection 

As shown in [TO] , in relation to the construction analogous to (|3.2p in the case of the Gaussian /3-ensemble, 
it is possible to deduce a three term recursion for the characteristic polynomial P n (x) and thus associate the 
inductive construction with a random tridiagonal matrix. The key idea for this purpose is to apply what 
may be regarded as the inverse operation of bordering, namely a random corank 1 projection 

T 

uu 

(4.1) IT„ + idiag A n+ iII„ + i, n n+ i = I n+1 - 77—^, 

\m\ 

for a suitable random vector u. 

The projected matrix will always have a zero eigenvalue. We want the remaining eigenvalues to be identical 
in distribution to diag A n . To determine the necessary form of u we note that we must have 



In+l — 



o n 

T 1 



Ah 



Ln+l 



o n 

T 1 



A n 

T 



Next we write 

A n+1 = Fdiagv4 n+ iV rt 

where V is the (n + 1) x (n + 1) unitary matrix of eigenvectors. We see thus immediately that we can choose 
u T = — iv T , where v T is the final row in V, and thus the vector of final components of the eigenvectors of 
A n +i- The structures (|3.ip and (|3.5p tell us that the latter can be chosen to have the form 

(iqi,-iqi, ■ • • , iq( n +i)/2, -*<7(n+i)/2), (n + 1) even 
(iqi,-iqi,...,iq n /2,-iq n /2,ic), (n + 1) odd 

where qi > 0, c > 0, normalized so that 

(n+l)/2 n/2 

(4-3) 2 £ ql = l, 2^+c 2 = l, 

i=i i=i 
depending on whether (n + 1) is even or odd. 



(4.2) 



Substituting — i times (|4.2[) for u T in (|4.2[) . the result of [TU1 Lemma 1] for the characteristic polynomials 
of the original matrix A n+ i, and the projected matrix, which by the choice f|4. 2|) has the same eigenvalues 
as A n , shows 



2 n/2 2 

'" 'ECT (n + l)odd 



(4.4) JnW_ = , 

[ ' AP n+1 (A) ^ 



A 2 j-* A 2 - A 

(n+l)/2~ „ o 



2_, \2 „ X 2 . + eVen 



2=1 



where {A^} are the positive eigenvalues of A n+ \. Rather than compute the distribution of {qf} and c 2 
directly, we will make a trial choice and show that it leads to the correct joint distribution of the zeros of 
P„(A) and P„ + i(A), which through (|4.4[) uniquely determines the distribution of {qf} and c 2 . 

Our trial choice is to have each 2qf distributed as 2wf in ()3.3|) . (|3.4|) . and to have c 2 distributed 
as b 2 in ()3.4|) . but with the further constraint of the normalization conditions (14. 3| . Since a normal- 
ized multivariate gamma distribution is a Dirichlet distribution D n [si, . . . ,s n ], our trial choice is that for 
(n + 1) odd {2 9 2 } i=1: ....„ /2 }U{c 2 } is distributed according to D n/2+1 [((3 /2) n / 2 , /3/4] (here {(3/2) n / 2 is short- 
hand for (3/2 repeated n/2 times), while for (n + 1) even {2g 2 } i=1 .... ! („ +1 )/2 is distributed according to 
£>(, i+1 )/ 2 [(/3/2)(™ +1 ' ) / 2 ]. In this circumstance, the conditional PDF of the positive zeros of P„(A) (and thus 
the positive eigenvalues of A n ) given the zeros of P„ + i(A), or equivalently the zeros of the random rational 
function in (]4.4|) . can be obtained by appealing to a known result, implicit in the work of Dixon [6] and 
Anderson [2], and given explicitly in [TQl sentence below (4.10)]. 

Proposition 4.1. Consider the random rational function 

n 

m = E 



1 A 



where each at is real and {ci} has the Dirichlet distribution _D n [si, . . . , s n ]. This function has exactly (n — 1) 
roots, each of which is real, and for given {ai} these roots have the PDF 



n — 1 n 

i=i j=i 



, , , \ n ( x i x j 

u kx r(siH + S n ) l<i<j<n-l -r-r T-r 

K 1 T( Sl )---T(s n ) II (ai-o^+^-i 1111 

l<.i<j<n 

supported on 

(4.6) ai > x\ > 02 > x 2 ■ ■ ■ > x„_i > a n . 

For (n + 1) odd, Proposition 14.11 applies to (|4.4|) with A i— * A 2 , Oi i— ► A 2 (i = 1, ... ,n/2), a n /2+i = 0, 
and Si = /?/2 (i = 1, . . . , n/2), s n / 2 +i = /3/4. For (n + 1) even, for application to (|4.4|) we require A i— > A 2 , 
aj I— > A 2 (i = 1, . . ., (n + l)/2) and {ci}i=i....,( n +i)/2 distributed as s, = (3/2 (i — 1, . . . , (n + l)/2). We can 
now read off the sought conditional PDFs. 

Proposition 4.2. Assume the validity of our trial choice of distribution for 2qf and c 2 . For (n + 1) orfrf, the 
PDF of the positive eigenvalues of A n , given that the positive eigenvalues of A n+ \ are \\, . . . \ n n, is equal 
to 

(47) 2"/»r((n +1)0/4) .f 2 - 1 n (*?-*?) ft 2 

1 • J (r(/3/2))(«-D/ 2 r(/?/4) 11 A 2 (3^/4-D ^1 1 ^ (A 2 - a 2 )^ 11 |x * ' 1 

w/iere 



(4.8) 



Ai > x\ > A 2 > • • ■ A„/ 2 > x n / 2 > 0. 
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For (n + 1) even, the PDF of the positive eigenvalues of A n , given that the positive eigenvalues of A n+ \ are 
Ai, . • • , A(„ + i)/2, is equal to 



(4.9) 



2("+ 1 )/ 2 r((n + l)/?/4) 

(r(/3/2))(» +1 )/ 2 



n 



(n-l)/2 11 ( X i ^j) (n-l)/2 (n+l)/2 

n * ^rr^§|3^ n n i*?-^- 1 



n 

l<*<j<("+l)/2 



w/iere 

(4.10) Ai > xi > A 2 > • • • A(„_i)/ 2 > ^(„+i)/2 > 0. 

Multiplying (|4.7[) and (|4. 8f> by ()1 .5|) (n n + 1) and (|1.6p (nwn + 1) respectively gives the joint PDFs 



n/2 n/2 

(4.11, a v/ j ' II (Af-ADCxf-xDn^-Air 2 " 1 , 



C„+i,/3 i=1 l<»<i<n/2 

7i even, subject to the interlacing (|4.8|) . and 



(n+l)/2 



^ n 



(n-l)/2 

n 



7/2 



(4.12) 



x n ^ - a 2 ) n (4 - 4) n w - 

l<i<j<(n+l)/2 l<i'<j'<(n-l)/2 i,J=l 



odd, subject to the interlacings (|4.10p . Here 



- (rW2)) (»-l)/2 r(/3/4) 

0^-0^ 2 ™/2r((n + l)/3/4) ! 



Cn+1,/3 — C n +i t f3 



(r(/3/2))(" +1 )/ 2 

2("+i)/2r((n + 1)^/4)' 



(n+1) odd 
(n + 1) even. 



On the other had the same joint PDFs can be obtained by interchanging the symbols Xj <-* Ai in (|3.7p and 
(|3.9p , and multiplying by (|1.5p (with Aj 1— » a;,-) and (|1.6p (with Aj 1— > Xj). This tells us that our trial choice 
of the distributions of the components of the eigenvectors is correct, and so in particular the qualification 
starting off the statement of Proposition 14.21 can be removed. 

The fact that our trial choice of the distributions of the components of the eigenvectors is correct also 
implies that we can make the replacements in (|4.4I) as indicated in the first sentence of the paragraph below 
(14.41). and so obtain 



(4.13) 



where 



(4.14) 



AP„+i(A) 



1/2 
^2+E 

i=l 
(n+l)/2 



A 2 _ A 2 > ( n + 1 ) odd 



E 



2wj 



A 2 - A 



2 • 



(n + 1) even 



1/2 



N n = < 



6 2 + ^2u; 2 , (n + 1) odd 

i=l 
(n+l)/2 

^ 2w 2 , (n + 1) even 



We can substitute (|4.13p with n i— > (n — 1) in (|3.3p , l|3.4p to deduce a random three term recurrence for 

{Pn(X)}- 
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Proposition 4.3. The characteristic polynomials {P ra (A)} n= 2 j 3,... for the matrices {A n } n= 2,3,... defined in- 
ductively in Proposition \S.1\ satisfy the random three term recurrence 

(4.15) P n+1 (\) = \P n (\)-b 2 n P n ^{\), n = l,2,... 

with -Po(A) '■= 1; Pi (A) := A and where b n has distribution T[nf3, 1]. This is the three term recurrence for the 
characteristic polynomial of the tridiagonal matrix 

Proof. The substitution gives 

Pn+lW _ i _ *T -Pn-l(A) 

AP„(A) AP„(A) 
and so we can obtain (I4.15|) with b 2 n — M n . The distribution of b 2 n then follows from (14.14p and the fact that 
the number of degrees of freedom in the sum of independent gamma distributed variables adds. That the 
three term recurrence relates to an anti-symmetric tridiagonal matrix is a standard result. □ 
Combining Proposition 14.31 with (|4.4p allows the distribution of the first element in each of the indepen- 
dent eigenvectors of (| 1 .4[) to be determined. We can thus show that the vector of first components of the 
independent eigenvectors of the random tridiagonal matrix (|1.4[) . which we choose to be positive, has the 
Dirichlct distribution given by Theorem 11.21 

First proof of Theorem 11.21 For a general n x n real symmetric matrix X we have that 



(4.16) 



Pn-l(A) _ Ci 



where P n (A) is the characteristic polynomial of X, P„_i(A) is the characteristic polynomial of the submatrix 
formed from X by blocking the first row and first column, {A.;} are the eigenvalues of X, and {c,} are the 
first component of the eigenvectors. For the matrix (|1.4p . the eigenvalues and eigenvectors have the special 
structures (|3.5p and (|3.ip respectively. Substituting in (|4.16p we reclaim again (|4.4p provided the independent 
members of {ci} are identified with the independent entries in (|4.2p (which are the last components of the 
eigenvectors of (|3.2p ). But we know from the paragraph below Proposition 14.21 that the latter entries have 
the Dirichlet distributions as claimed. □ 

5. Method II: Jacobians of anti-symmetric tridiagonal matrices 

In this section, we present an alternative proof of Theorems 11.11 and 11.21 based on the mapping between 
a real anti-symmetric tridiagonal matrix and its positive eigenvalues and first row of the eigenvector matrix. 
This proof is very much in the spirit of [7]; in fact, many of the results used there for symmetric matrices can 
be used here for anti-symmetric ones, with very minor modifications. In the interest of brevity, we present 
the proofs only if the modifications are non-trivial. 

We start by giving a set of linear algebra results, which build up to the computation of the eigenvalue 
PDF for the random tridiagonal matrices (|1.4|) . 

Anti-symmetric matrices are normal matrices, i.e., they have the property that they commute with their 
transpose (for a treatment of normal matrices, see any linear algebra book, e.g. |13j). Equivalently, they 
have the very special property of being diagonalizable via a unitary transformation. Any real anti-symmetric 
matrix T can be decomposed as T — UAU H , with U a unitary matrix and A the diagonal matrix of 
eigenvalues. 

Let T be a real anti-symmetric tridiagonal matrix in reduced form, defined as 

&„_i 
-b n -i 6„_ 2 



(5.1) T 



6i 
-h 



if bi > for all i — 1, . . . , n, then when n is even T is full-rank, whereas when n is odd is a simple eigenvalue. 

As mentioned already in the introduction, the non-zero eigenvalues for such matrices come in pairs 
(i\j, —i\j), with j = 1, . . . , , and we assume the ordering Ai > A2 > . . . > Aj„j > 0. 
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The decomposition T = UAU H is not unique; to make it unique, we impose two conditions. First, we 
order the diagonal of A as follows: ■ ■ ■ , iAr«i , — iAi, . . . , — «Ar„i ). If n is odd, we let A(n, n) — 0. (Note 

that this ordering is different to that used in Proposition 13. 11 but is more convenient for present purposes.) 

Second, we impose the condition that the first row of U has positive entries (note that eigenvectors are 
only defined up to multiplication by rotations e lB , and also that if Uj is an eigenvector for eigenvalue iXj, 
then 7ZJ, the conjugate of Uj, is an eigenvector for — iX). 

The first result we state is very similar to Theorem 7.2.1 in [T7], and the proof given in [T7] is adaptable 
almost verbatim to this case (thus we choose not to repeat it). 

Lemma 5.1. With the notation above, let q := (gi, . . . ,q n ) be the first row of the matrix U , and X be the 
diagonal of — iA. Then T is uniquely determined by X and q. 

Due to the real anti-symmetric nature of T, we can in fact deduce that 

Corollary 5.2. T is determined by Ai, . . . , Ajnj , qi , . . . , <?j » j . 

Proof. Note that by our choice of U and A, gj„ = qj for all j = 1, . . . , . The statement is immediate 
for n even; if n is odd, note that is also an eigenvalue and that the first component q n corresponding to it 

n 

can be determined from the condition q 2 = 1. □ 

3=1 

Before we proceed, we will need the following lemma, which is similar to Lemma 2.7 of [7]. 

Lemma 5.3. Let A(A 2 ) := JJ ( iA; — X 2 ), where X\, . . . , Aj„j are the positive imaginary parts of the eigen- 
values of T , and q is as above. Then 
• if n = 2k. 



n-l 

n H 



(A(A 2 )) 2 



2* n fe 2 Ao 

1=1 



• ifn = 2k+ I, 



(A(A 2 )) 5 



n-l 

n % 

i=l 



n ^ 



Proof. To keep eigenvalues real, we will examine the matrix iT. Just as in [7], we use the three-term 
recurrence for the characteristic polynomial of iT. Denote by X\ m \ i — 1, . . . , m the eigenvalues of the 

m . . 

m x m lower corner submatrix of iT, and denote by P m (x) — Y[ ( x ~ X ™ 1 ) the corresponding characteristic 

i=l 

polynomial. Then the three term recurrence (|4.15[) holds, and from this, writing 

m rn — 1 



(5.2) 

we deduce that 
(5.3) 



n 



A„ (rn) - A 



(m-l), _ 



1<><K 
l<j<m 



ni^ m -i(ADH n i p ™( A 



(m-l) | 

3 



m — 1 



n Pm( - x 



(m-l). 



= b 



2(m-l) 



n p ™- 2 ( a 



(m-l) 



By repeatedly applying (|4.15p and (|5.2p . we obtain that 



n-l 



(5.4) 
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Like in [7], we make use of a simple identity for qf, which is a particular form of Theorem 7.9.2 from [17) : 



(5.5) q\ = P "~ l( ^" ) , V 1 < i < n . 

Note that (|5.5[) can also be seen as a corrolary of (|4.4[) . 
Let us now examine P^(Xi), with A^ = \\ ■ Since 

n 

»=l jV» 

it follows that 

k 

• if n = 2k, |i^(Ai)| = 2|Ai | J] |Af - Af j, for alH = 1, . . . , k, and 

3=1 

• ifn = 2fc + l, 



o 



P^(Xi)\ = 2Af II | A? - Af |, for alH = 1, . . . , k, and 



3 = 1 



A- 



° i^(o)i = n a 



2 



i=l 



Since = Qj? f° r all j = 1 , . . . , [^] , it follows that 

• if ra = 2fc, 



niPn-i(Ai)i n&? 

i=l i=l 



2«(A(A2))4 ft Af 2«(A(A2))4 ft Af 



(5.6) 

• if n = 2k + 1, 

n n 

n i^-^aoi n ^ 
(5.7) *MU 4 - i=1 



i=i 



2«-i(A(A 2 )) 4 ft Af 2«(A(A 2 )) 4 ft Af 



i=l i=l 

Rewriting (|5.6[) and (|5.7[) and taking square roots, one obtains the statement of the lemma. □ 

Consider now the transformation T <-» (g>, A), with all of the conditions imposed above; the transformation 
is one-to-one and onto, and it must thus have a Jacobian J . We will compute this Jacobian in the same way 
we computed the Jacobian of the similar transformation for symmetric tridiagonal matrices in [7j: we will 
make use of the fact established in Section [3] that the anti-symmetric GUE ensemble has the same eigenvalue 
distribution as the tridiagonal model (|1.4[) with (3 = 2. 

An alternative derivation for the Jacobian J is given in the Appendix. 

Specifically, we will require the following three properties of the anti-symmetric GUE ensemble (J3 = 2). 

Property 1. The eigenvalue distribution of the anti-symmetric GUE ensemble is given by (jl.4[) with (3 = 2. 
Property 2. The set of eigenvalues and the set of eigenvectors of the anti-symmetric GUE ensemble are statistically 
independent of each other. 

Property 3. • For n = 2k, let U = [ui, U2, . . . , u^., u7, uj, . . . ,Tife] be a matrix of unit-norm eigenvectors. For 
each j = 1, . . . , k, write Uj = Vj + iv/j, with v, and real. Then, with probability 1, the set 
{Vj-, j = 1, . . . ,k} U {Wj, j = 1, . . . , k} is an orthogonal basis for 0(2k). Moreover, for any j, \/2vj 
and -\/2wj are distributed uniformly over the sphere. 

• For n = 2k + 1, let U = [ui,U2, .. . , u^, TTf, U2, ■ ■ ■ , Ufc, z] be a matrix of unit-norm eigenvectors. 
For each j = 1, . . . , k, write Uj = Vj + i~Wj, with Vj and Wj real. Then, with probability 1, the set 
{vj, j = 1, . . . , k} U {w^, j = 1, . . . , k} U {z} is an orthogonal basis for 0(2k + 1). Moreover, for any 
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j, v2~Vj and \^2wj are distributed uniformly over the sphere, and z is distributed uniformly over 
the sphere. 

Choose now U to be the (unique!) matrix of anti-symmetric GUE eigenvectors which has on its first row 
all positive numbers q\, . . . ,q n . From the three properties above one can immediately deduce Proposition 



Proposition 5.4. If n — 2k, (qi,...,q n ) has the same distribution as rr^^jrr , where w is a vector of k 
independent variables with distribution w ~ (\2, X2, ■ ■ ■ , Xi)- 

If n = 2k + 1, (qi, . . . , q n ) has the same distribution as iT^^^ru > where w is a vector of k independent 
variables with distribution w ~ {xii Xii • ■ ■ > X2), and z is a scalar Xi~ distributed variable independent of w. 

We now proceed to compute the Jacobian of the transformation T <-* (q, A). 
The joint element density of a matrix T from the distribution (|1.4p with (3 = 2 is 

n-i -"fl b2 



(5-8) Kb) = — Il^e 



i = l 



nr(f) 



i=l 

,2/ 

i=l 



Denote by db = A^L/d&i, dA = A^jAj. To be consistent with Property 3, denote by dq 

• if n = 2k, dq is the surface element on the fc-dimensional sphere of radius 1/2, 

• if n = 2k + 1, dq is the surface element on the (k + l)-dimensional ellipsoid of first k axes equal to 
1/2, and the (k + l)st equal to 1. 

With the transformation T (q, A), it follows that 

(5.9) K b ) db = J KKQi a )) d <? dA = HQ, A ) dgdA . 

By Proposition [5TU and Properties 1, 2, and 3, it follows that 

• if n = 2k, 

k 



(5.10) u{q,X) = -L_ (A(A 2 )) 2 e~ M JJ* i 

• if n = 2fc + 1, 



(5.11) u{q,X) = j^— (A(A 2 )) 2 HA? e"« 1 -^TT 1 II » ■ 

C ^ 2fe + 1 f=l r (2) f=l 

Hence all that is left is to compute the Jacobian J as 

= Kg, A) 
Kb) • 

Since the Frobenius norm of a matrix is preserved by orthogonal similarity transformations, 



(5.12) f> 2 = i||T|| F = i||A|| F = ^A 



^2 

i=l 



By putting together (|5.8[) - (|5.12p . Lemma T5.3[ and the definition of Ci. n from Theorem II. 1[ all constants 
cancel and we obtain the following lemma. 

Lemma 5.5. The Jacobian J of the transformation T <-> (q, A) is given by 
• if n = 2k, 

n-i 

n bi 



n 

1=1 
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ifn = 2k+ 1, 



J 



n-l 

n f>i 

i=l 
k 

i=l 



We are now able to give an alternate proof of Theorems 11.11 and 11.21 
Second proof of Theorems \1.1\ and \l.^ Starting from the joint element density of the matrix A^, 



(J>n,p(b) 



nr(f) 

i=l v 7 



¥-1 

e 1=1 



we make the transformation «-> (g, A) , into the eigenvalues and first row of the eigenvector matrix; we 
use the computed Jacobian from Lemma 15.51 and the expression for the Vandermonde from Lemma [ 



For n = 2k, we obtain the joint density of q and A: 



Vn,0{q, A) 



JVn,/3(b(q,\)) 

2 n-l 



n — 1 ip 

— e ,=1 



nr(|) ru 

i=i v 7 »=i 



A, 



i=l 



n — 1 / \ 

nrf 

i=l v 7 
2n-l+fe/3/2 
n— 1 / \ 

nrf 



2 fe n 



0/2 , k N /3/2 

2 fe n ?|Ai 



1=1 



2^ TT i/ 3 / 2 - 1 



(A(A 2 ))^ n A « 



fe 

n 

i=l 



- E A? 
e i=1 



V.i=l 



|8-1 



It is easy to check that 



2n-l+fc/3/2 -y 



2fe-l+fe/3/2p 



(f) 



n'r(f) c »'» (r(f 



this shows that A and q decouple and that they have the distributions described in Theorem 1 1.1 1 and 

o 

is 



For n = 2k + 1, the joint density on q and A can be obtained as 
v n,p{q,\) = JVn,p(b(q,\)) 



n— 1 ip 

n v 



n-i / v fe 

nr | qnUQi 

1=1 v 7 1=1 



- E A? 

e i=1 



' n &; * 

i=l 



n— 1 



nr(f 

1=1 v 7 

2n-l+fc/3/2 

= n'r(f) 

i=l v 7 

Once can easily check that 



2 fe g„ n 8? A 3 



/ 



/3/2 / fe N 0/2 

2 fe 9« II 9?A? 

t=l 
fe 

<7n I! 



E A? 



i=l 

(A(A 2 )) ? nAf/» e "£ A| 



^ 2 fe+fc/5/2 p ^ fc/3+/3/2 ^j 



n—1 / \ 

nr(f 

1=1 v 7 



/3,2fe+l F //3 



r(f ) (r 



(«) 



so again, A and g decouple, and they have the distributions described in Theorems 11.11 and 11.21 
This completes the proof. □ 

6. Method III: an orthogonal transformation 

This proof is based on the observation that the distributions of Theorem 1 1 . 1 1 are . essentially, /3-Laguerre 
distributions: they are the same as the distributions of singular values of the -B^n.a, bidiagonal, root-Laguerre 
matrices given in [7]. 

We describe this class of matrices below. 

For 2a - (n- l)/3 > 0, let Bp. Ui a be the n x n random bidiagonal matrix 



B/3,n,a 



X2a 

X(n-l)f) X2a-P 



XP X2a-(n-l)P 

It was proved in 7] that the eigenvalue distribution of Lp jTlta — Bg n a Bg n a (or, equivalently, Wp, n ,a = 
.Bj n a B g n a , since they are the squares of the singular values of Bp^ n<a ) is the well-known Laguerre distri- 
bution of size n and parameter a — (n — l)/?/2 — 1. Denoting by Ai > . . . > A„ > the eigenvalues of Lg iTlta , 
their joint PDF is given by 



o-(n-l)|-l "EAi/2 

2 e i=1 



where 



2 -na j~J 



(«) 



Since the singular values of Bp. rha are the square roots of the eigenvalues of Lp in<a , the joint PDF of the 
singular values of Bp tTlta is given by 



(6.1) 



fn ,a,P 



i<j 



10 



The proof is based on the observation that, if we let the size of the matrix anti-symmetric matrix be 
n = 2k, with k variables, the parameter a — 2 7 1 (3, and in addition scale each <Ji by \[2, the PDF of (|6.1[) 



is the same as the one in Theorem ll.il while, if we let n = 2k + 1 and 
the PDF of (|6.ip is once again the same as in Theorem 11.11 

We first recall the following well-known result in linear algebra. 



_ 2fc+l 



P, after scaling Cj by v2, 



Proposition 6.1. Let Y be a real n x n matrix, and construct the 2n x 2n matrix 



V Y 







-Y 




Let o~i, . . . ,o~ n be the singular values ofY (with multiplicities) . Then the eigenvalues ofVy are ±io~i, ±io~2, ■ ■ ■ , ±io~„ 
(also with multiplicities). 

If Y is a bidiagonal matrix, one can "shuffle" the entries of the matrix Vy to make an (anti-symmetric) 
tridiagonal out of them. We first need to define "shuffling". 

Definition 6.2. We define the "perfect shuffle" 2n x 2n permutation matrix P n to be given by 

1, ifj = ^r orj = n+i V2n>i,j>l 



0, otherwise. 



Note that, given a matrix X, P n X has the same rows as X, but listed in the following order: 1, n + 1, 2, n + 
2, . . . , n, 2n, while XP^ has the same columns of X but rearranged in the same order 1, n+1, 2, n+2, . . . , n, 2n. 
Also note that, since P n is a permutation, P n is orthogonal. 

We can now define the alternating sign perfect shuffle (ASPS) matrix. 

Definition 6.3. Let D n be the diagonal matrix for which D n (i,i) = (— 1)W, i = 1..2n. We call the matrix 
Q n — D n P n the alternating sign perfect shuffle (ASPS) matrix. Note that Q n is also orthogonal. 

We can now explain the effect of the ASPS matrix on a tridiagonal matrix. 

Lemma 6.4. Let T be a tridiagonal anti- symmetric matrix (as in (|5.1|) ). Then T — QnVsQn, where B is 
the n x n bidiagonal matrix 



B = 



K 

b n -i b n -2 

b n -3 



b n - 



bi 



and Vb is like in Provosition \6 . 1\ 



The proof of Lemma 16.41 is an easy exercise; it suffices to see how the entries of Vb move around under 
left multiplication by Q n , respectively, right multiplication by Q^. 

We give here an example: for i < n, the entry (i, i + n) moves first to (2i — 1, i + n) under the multiplication 
by Q n to the left, then it moves to (2i — 1, 2(i — 1)) under multiplication by Q„ to the right. Along the way, 
it gets multiplied by (— ^M + I"] , and thus it changes sign. 

The other cases can be examined in the same way. 

The proof for matrix (A^) size n = 2k differs slightly from the one for n = 2k + 1; we present them 
separately. 

Third proof of Theorems 11.11 and 11.21 n — 2k. Armed with Lemma f6T4| Theorem 11.11 follows directly in 
the case when n — 2k, as 

\f2Ai = Q n V B ^ a Ql , 

where the equality should be understood in terms of distributions. 

In addition, the distribution of the first components of the eigenvectors of A^ (given by Theorem 1 1.2[) can 
be obtained from this orthogonal similarity transformation, as a consequence of the following three facts: 
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if B = UYV T is the SVD of B, then 



V B 



U 
V 



u 

-V 



-E 
E 



U T 



V T 
-V T 



is the eigenvalue decomposition for Vb- 

• the first row of the left singular vectors for Bp i7ltCt is distributed like a (normalized to 1) vector of 
i.i.d. xp random variables (see [7]); 

• the first row of the matrix Q n X is the same as for X, for any matrix X . 

□ 

For the n odd case, the proof is only slightly more complicated. 
Third proof of Theorems 11.11 and II .21 n = 2k + 1. The first obstacle in using the ASPS matrix is that 
the size of A@ is odd. This is easily overcome by introducing an extra row and column of zeroes, set 



A 2k+l °2fe 







where 02k is the column vector of 2k zeroes. We immediately obtain that 
where 

Xkf3 

Xih^ip X(k-t)p 

X(k-2)/3 



C/3,fc = 



X.E 







with equality here being in the sense of distributions. 

This is a second obstacle, as Cp^ is not a /3-Laguerre matrix, and has as a singular value. The latter 
part can easily be corrected by removing the last column of Cp t k and creating a (fc + 1) x k matrix Cp } k- 

We would now like to show that the k x k matrix Lp ^ — Cj k C/3^ has the same eigenvalue distribution 
as the matrix Wp.k, a = B^ k a B p k a with a = ^f^p. 

Notation. We will now make the following notational convention: in the below, any variable indexed by i 
(e.g. a,i) will have distribution Xi/3/2- Some indices will therefore be skipped. 

If we denote the entries of Cp.k as follows: 



C/3,fc 



b2k 

&2fc-l 



?2(fc-2) 



then 



b 2k + b 2k-l hk-lhk-2 

&2fc-l&2fc-2 Hk-2 + Hk-3 &2fe-4&2fe-3 

b 2 b 3 b\ + b\ 



while at the same time, if we denote the entries of Bp^,a by 



B 



f3,k,a 



&2k-2 d2k-l 



0,2 a 3 
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and 



l 2k+l 



a. 



2k-2 



a2fc-lfl2fc-2 



a2fe-l02fc-2 
2 2 
°2fc-l + a 2fc-4 2 fc-4a2fc-3 



0203 



Note that, with the notational convention adopted above, the marginals of the entries of Lp t k an d Wp_k,a 
are the same, since independent chi-square variables add to a chi-square variable. 

Claim 6.5. The Choleski factorization of the matrix Lp,k yields a matrix whose distribution is the same as 

Wp,k,a- 

Proof. Note that Cp^ is not the Choleski factor of Lp f., because it is (k + 1) x k instead of k x k. 
We will prove that if we solve the system of equations 



c 2i+l + ^21-2 
x 2i+l x 2i 



b^ + b^-i , fori = 1,2,. 



hi+ib2i , for i = 1,2, 



. . , k , 
k-1 



with (bf, 



b 2 

' u 2k+l 
[X 2 , . . . , x 2k-l) ~ (^2j 



being independent chi-squared variables of parameter ((3/2, . . . , (2k + l)/3/2), then 
Moreover, we will obtain as a bonus that aj|fe+i ^ s chi-square distributed, 
independently of all others, with parameter 2 /?. 
We first need the well-known lemma below. 

Lemma 6.6. If x ~ Xr an d V ~ Xs; anc ^ x i 2/ are independent, then z — is distributed like Beta(r, s), and 
z is independent of (x + y). Moreover, if w ~ Xr+s independently of x and y, then wz ~ x 2 , w(l — z) ~ 
and is independent of w(X — z). 



First, we find a; 2 , and £§: 



&1 + &2 , 

6 2^i_ 

3 b\ + b 2 



It follows immediately from Lemma 16761 that (63 — x 2 ,, x^x 2 ,) 



2 

~ X 2i-2> 


we can obtain 


T 2 


= & + 






X 2i 






^21+1 



t '2j-2 



and 



Assume now that 



(b 



2 

2i-l 



2 2 2 

x 2i-2i x 2i x ?>i 



i J '2i-1 



for some i > 2. From the formulae above and Lemma \6. 61 one sees that (x 2 ,^ x\ i+l ) ~ (6^, bfi+i) ano - that 
they are independent of all x 2 with j <2i — 1. Furthermore, is independent of all a; 2 - with j < 2i + 1. 

Altogether, this yields 



(^ 2 



_ 2 22 

2i+l x 2i > a; 2 j x 3 i • 



' 1 ■ h 2i+\i 



2 »,2 l2 
2j °3) • 



one can easily see that x?, k+1 



b\k 



b\ and that it is independent of all other b's, and thus of all other x's. 



Thus the claim is proved by induction. 



□ 



It remains to conclude that, since the Choleski factorization of the matrix Tp,k yields a matrix whose 
distribution is the same as Wp^,a, Theorems 11.11 and 11.21 are true for n = 2k + 1. □ 
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7. Sturm sequences and Prufer phases 



The characteristic polynomial P n (X) for %T with T the anti-symmetric tridiagonal matrix (|5.1[) is identical 
to the characteristic polynomial for the symmetric tridiagonal matrix T s obtained from T by removing the 
minus signs below the diagonal. This can be seen from the three term recurrence (|4.15l) . Hence some 
fundamental results applying to the characteristic polynomials of real symmetric tridiagonal matrices also 
apply to P n (\). 

One such result relates to A(/i), the number of eigenvalues less than \i. With := — Pi(/i)/Pi_i(/i) 
(i = 1, . . . ,n) the theory of Sturm sequences (see e.g. E]) tells us that N(/j) is equal to the number 
of negative values in the sequence {ri}i=x <n (assumming that fj, is not an eigenvalue of T). In our anti- 
symmetric setting, there are precisely [{n + 1)/2] eigenvalues less than or equal to 0; since [(n+ l)/2] is equal 
to the number of terms in {r<}i=i 3,..., giving the following result. 

Lemma 7.1. Let A r+ (/i) be the number of positive eigenvalues of iT less than or equal to fx, and suppose /i 
is not an eigenvalue. We have that N + (fi) is equal to the number of negative values of {?"j}i=2,4,... minus the 
number of positive values of {ri}i=i,3,.... 

This result relates in turn to shooting eigenvectors x = (x n , . . . , x\) T of the symmetric tridiagonal matrix 
T s . With xi and /1 given, the shooting eigenvector is specified as the solution of all but the first of the n 
linear equations implied by the matrix equation (T s — fj£)x = 0. Further, with x n +i defined as the first 
component of (T s — /il)x, we have that Xi/xi-i = — ri_i/&i_i as can be checked from the recurrence (|4.15|) . 

Finally, we discuss the Prufer phases 0f associated with the shooting vectors of T s (see e.g. [T5]). For 
2 < i < n these are specified in terms of the characteristic polynomial by 

/» i\ . a u 1 P>-i(m) 

(7 - 1} COt "' Cr^M 

where b n := 1, together with the condition that 

(7.2) 0£ -*• -[?/% as 

(this is consistent with (17. 1|) as the RHS — > 00 when /i — > 00) and the requirement that (9f be continuous in 
Differentiating (|7.ip with respect to fj, gives 



Vd M y s in 2 ^ 6?_ X V (^- 2 (M)) 2 

But according to the Christoffcl-Darboux summation formula the RHS is positive so we recover the well 
know result (see e.g. [TS]) that 0j (j > 2) is a strictly decreasing function of /1. Further, we see from (|7.ip 
and (|7.2[) that 6*f = 7r/2 + fc7r, k = 1, 2, . . . for the fc-th positive zero of Pj_i(^t), while 6*f = /c7r, /c = 1,2,... 
for the fc-th zero of Pj_2(/i). Since ^ = is a zero of Pj(fi) for j odd, it follows in particular that 0° ■ = 7r/2, 

fl^-i - 0. 

In a significant advancement [Jj5] (see also [H>]), the Prufer phases associated with the tridiagonal matrix 
(jl.ip . p.2p relating to the Gaussian /3-ensemble, have been shown to satisfy a stochastic differential equa- 
tion. An analogous study of the Prufer phases of the present anti-symmetric Gaussian /3-ensemble awaits 
investigation. 
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Appendix 

As mentioned below (|5.7p , the strategy of our proof of Lemma 15.51 was first used in [7] to compute the 
Jacobian for the change of variables from the elements of a general real symmetric matrix to its eigenvalues 
and first component of its eigenvectors. Subsequently, a direct approach to the computation of this Jacobian 
was given In this appendix, we will show how this direct approach can be adapted to provide an 

alternative proof for Lemma 15.51 

Suppose first that n is even. The starting point is the identity (|4.4[) as it applies to (|5.ip . with the LHS 
rewritten to read according to Cramer's rule. Thus, after further minor manipulation, we have 

n/2 2 , 



(Al) ( (I "- Air) ~ 1 )n = ?l 



3=1 3 

By equating successive powers of A on both sides we deduce 

n/2 n/2 n/2 

1 = £ 2ql bl_ x = H X l * + 6 n-2^n-l = E 

3=1 3=1 J=l 

n/2 n -l n/2 

(A.2) * + taW-. = E 2 <^, • • ' • * + II b " = E 2 «l A f ~ 2 

3=1 i=l 3=1 

where the * denotes terms which have already appeared on the LHS of preceding equations. In particular, 
the set of equations (|A.2|) is triangular in 6„_i, 6„_2, ■ ■ ■ ,b±, and the first of these equations implies 

n/2 

(A. 3) q n /2dq n /2 = - E q 3 dq J ■ 

3=1 

Taking differentials of the remaining equations, substituting for q n /2dq n /2i an d taking the wedge product of 
both sides shows 

n-l 9 n-l "/2 

( a - 4 ) n $- x * = — n det - ^ bA^- 1 ] 

i=i 3=1 

Here, to obtain the LHS essential use has been made of the triangular structure, and we have written 

n/2 n/2 

(A.5) db := h^db 3 , dX := f[ dX 3 , dq := J] d qj . 

3 = 1 i=i 

According to [TTJ Proposition 2.1], up to a sign the determinant in (|A.4[) evaluates to 

m/2 

n (^(A 2 

3=1 

so after making use too of Lemma 15.31 we have 

„ i n™=i t>j 

db = — dX A da. 

^n/2Y[%\q 3 X 3 

Noting that dq L /2q n / 2 — dq, as follows from the meaning of dg below (|5.4|) and a simple scaling, we read off 
the Jacobian J as specified in the first case of Lemma 15.51 
For n odd, rewriting (|4.4p as in (|A.1|) gives 

(«-l)/2 9 2 

((I„-AiT)-) n = E H^+ C 2 



, 4 

\2\ 



3=1 J 
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We see that the first equation in (|A.2[) needs to be modified to read 

(n-l)/2 

1= E 2 ^ 2 + c2 > 
j=i 

and the remaining equations remain the same except that the upper terminal in the summation on the RHS 
must have n/2 replaced by {n — l)/2. In particular, the variable c does not appear in any other equation, 
and we obtain in place of (|A.4[) 

n-1 ("-l)/2 

(A.6) H bf -'db = 2"- 1 Yl gfdet 

j'=i i=i 

where dA and cZq are as in (|A.4[) but with the upper terminals in the products replaced by (n — l)/2. The 
argument used to establish 11, Proposition 2.1] shows that up to a sign the determinant is equal to 

™/2 4 

ii a K a(a2) ) • 

3=1 

Hence, after substituting in (|A.6p and using Lemma l5.3l we obtain for the Jacobian 

c 2 nfe 1)/2 *;V 

This agrees with the second case of Lemma l5.5l after noting that dq/c = dq, with dq as specified below (|5.4|) . 
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